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ABSTRACT 

We present first results from electronic Multi- Element Remotely Linked Interferometer 
Network (e-MERLIN) and electronic European VLBI Network (e-EVN) observations 
of a small sample of ultra-steep spectrum radio sources, defined as those sources with 
a spectral index a < —1.4 between 74 MHz and 325 MHz, which are unresolved on arc- 
second scales. Such sources are currently poorly understood and a number of theories 
as to their origin have been proposed in the literature. The new observations described 
here have resulted in the first detection of two of these sources at milliarcsecond scales 
and show that a significant fraction of ultra-steep spectrum sources may have compact 
structures which can only be studied at the high resolution available with very long 
baseline interferometry (VLBI). 
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1 INTRODUCTION 



Radio sources with an ultra steep radio spectrum (USS; 
defined here as having a < —1.4 where oc v") are of- 
ten associated with extreme phenomena in the universe. 
USS sources with a classical FRIT type radio morphology 
(highly coUimated jets and well-defined lobes with promi- 
nent hotspots) are, for the most pa rt, associated with very 
distant galaxies (2 < z < 5.2; IPe Breuck et all l200d : 
iRottgering et all Il997l ) . Being very massive and often lo- 
cated in proto-clusters, these sources can be used to study 
the origin and evolution of massive galaxies and clusters. 
Recent high-redshift results include (i) the finding that mas- 
sive radio gal axies with masses of around 10^^ Mq exist at 
up to 2 ~ 5 (jSevmour et al.l 120071 ) . and (ii) the discovery 
of proto-clusters up to z ~ 4 with total masses as ex- 
pected for progenitor s of nearby clusters (a few 10^* M©; 
IVenemans et al.ll2007l ). USS sources with mostly diffuse ra- 
dio emission are virt ually always associated with close-by 
clusters up to z ~ 0.5 l|van Weeren et al]|2009l 'l and the emit- 
ting regions can have spectacular sizes up of to a few Mpc. 
The most likely explanation for this class of object is ra- 
dio emitting plasma tracing shocks in merging clusters. The 
study of these two types of USS source started in the late 
seventies with studies of their morphology, and follow-up 
work with large optical and X-ray telescopes subsequently 



gave an understanding of their enormous importance for 
studies of galaxy and cluster evolution. 

A third class of USS has not yet received much atten- 
tion; those which are unresolved on arcsecond scales. The 
aim of the observations presented here is to take a first step 
towards determining the nature of these objects. There are 
a great number of possible scenarios for the origin of such 
compact USS objects including (i) radio galaxies located at 
or near the epoch of re-ionisation (it has been suggested 
that objects with the steepest sp ectral index correspo nd to 
the highest redshift objects, e.g. iKrolik fc Ch^ll99ll ). (ii) 
young obscured radio galaxies with the early phase perhaps 
accompanied by a starburst (|Reuland et al.l |2004| inferred 
star formation rates of 1000 Mq /year), (iii) steep spectrum 
core quasars with the steepening perhaps caused by cooling 
of the plasma via scattering of infrared photons from the 
torus by the synchrotro n-emitting electrons, and (iv) pulsars 
(IDe Breuck et al. I I2OOOI found a clear over density of sources 
with a < —1.6 near the Galactic plane). Long baseline radio 
interferometric observations may give very important clues 
about the nature of these objects. Star formation activity in 
galaxies is expected to p roduce low brightness te mperature 
emission (~ 10" K, see lAlexandroff et all |2012| and refer- 
ences therein) and diffuse morphology. Synchrotron emission 
from radio lobes on kiloparsec scales can be very well imaged 
with e-MERLIN, but will still be mostly resolved out on mil- 
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liarcsecond scales. Finally, steep spectrum quasars and pul- 
sars will be compact when probed with VLBI; the potential 
Galactic origin of the source can then be verified by multi- 
epoch proper motion measurements. This letter describes 
e-MERLlN and e-EVN observations of a sample of compact 
ultra-steep spectrum sources, aimed at attempting to inves- 
tigate their morphologies and narrow down the possibilities. 
Section [2] describes the sample selection and observations, 
the results are discussed in section [3] and in section |4] we 
draw some conclusions. 



2 SAMPLE SELECTION AND OBSERVATIONS 

The initial sample of USS so urces was found by correlating 
the 74 MHz VLSS catalogue ('Cohen et al."2007) against the 
325 MHz WENSS database ( Rcngclink ct al. 1997) and se- 
lecting the small fraction of VLSS sources with a sufficiently 
steep spectrum. A sub-sample of compact sources was then 
selected by choosin g objects which we re both unresolved in 
the FIRST survey jBeckor et al."l995'), and had a flux ratio 
between the NVSS |Condon ct aL 1998.) and FIRST cata- 
logues of 0.8 < R < 1.2, ensuring that the 1.4 GHz emission 
is indeed compact on scales of liarcsecond. The FIRST 
and NVSS maps were also visually inspected to ensure that 
the sources were unresolved and that there were no nearby 
objects which could have led to ambiguities in the measured 
fluxes. The resulting sample contained 17 sources. A sub- 
sample of five sources was selected for an exploratory VLBI 
survey, all of these sources had NVSS flux densities larger 
than 10 mjy (see Table [1} . 

2.1 e-EVN and WSRT observations 

The five sources (J072212-h291042, J082916+383453, 
J130612-H514407, J131655-h483200 and J151229-H470245) 
were observed with the EVN at 1.6 GHz in e-VLBI mode 
Hipmoru 2008) over two 10-hour runs carried out on 
10*'' and ll"' June 2010 (programmes EP070A and B; PI 
Rottgering). Ten stations participated in the experiment: Ef- 
felsberg, Medicina, Onsala, Torun, Westerbork (12 antennas 
tied array), Lovell, Cambridge, Darnhall, Knockin and She- 
shan. The target sources were phase-referenced to nearby 
calibrators; the typical on-source time was 2.5 hours. The 
total aggregate bitrate per telescope was 1024 Mb/s, except 
for the MERLIN antennas which participated with a lower 
rate of 128 Mbps. Both LCP and RCP polarizations were 
observed with 2-bit sampling. The data were processed us- 
ing standard procedures following the EVN Data Analysis 
Guide{3- 

During the e-EVN observations, synthesis array data 
were also recorded at Westerbork and were independently 
processed in AIPS. The data were recorded using 8 bands 
of 20 MHz, each split into 64 channels. For amplitude cal- 
ibration we used 3C286. There was significant interference 
during the observation resulting in the loss of the first two 
subbands. 

Out of the five objects from the sample which were ob- 
served with the e-EVN, J072212+291042, J130612-f 514407 

^ http: //www. evlbi.org/user_guide/guide / userguide.html 



and J151229+470245 were detected. Note that initially 
J130612-I-514407 was not detected because of a 12.8- 
arcsecond error in the pointing position; a clear detection 
of the source was obtained only after this error was discov- 
ered and a wide-field image was made. Table [2] summarises 
the detections of these sources. Since, for the sources which 
were detected, the e-EVN observations recovered only ~20% 
of the VLA fiux at the same frequency, e-MERLIN commis- 
sioning observations were requested in order to investigate 
their structures on intermediate scales. Comments on each 
of these sources are given in the next section. 

2.2 e-MERLIN observations 

e-MERLIN observations of two of the EVN-detected sources 
were requested in order to investigate their structures on in- 
termediate scales. The e-MERLIN observations were carried 
out in March 2011, during the commissioning period. At this 
point in the commissioning phase, five stations of the array 
were available with the new C-band (4- to 6-GHz) system: 
Mk2, Tabley, Darnhall, Knockin and Defford. Observations 
were carried out at 6.6 GHz using the new capabilities of the 
upgraded e-MERLIN arrajjf] and used four adjacent bands of 
128 MHz, each split into 512 channels, starting at 6.64 GHz 
with a central frequency of 6.89 GHz. In both cases 3C286 
was used as the primary flux calibrator and OQ208 was used 
as the bandpass calibrator. Both sources were observed for 
~4.5 hours on 26th March and ~10 hours on March 30th 
2011. J151229-f470245 used the source 1500-^478 for phase 
referencing, J072212-I-291042 used 0736+299 as the phase 
calibrator. Of the total observing time, only ~9 hours on 
each source were ultimately usable due to a variety of com- 
missioning issues. The data were flagged spectrally, averaged 
in frequency, fringe fltted and then calibrated using standard 
methods for phase referencing experiments. 



3 RESULTS 

3.1 J072212+291042 

With a flux density of lOmJy in FIRST and 11 mJy in the 
NVSS (TablelHl, J072212+291042 is the weakest of the three 
VLBI detections. The e-EVN observations showed a ~2 mJy 
source with no clear evidence of structure on VLBI scales, 
consistent with a point source model, albeit at signiflcantly 
lower flux density than that detected by the VLA and West- 
erbork: only ~30% of the Westerbork flux is recovered by 
the e-EVN at the same frequency indicating the presence of 
significant structure on intermediate scales. 

The source is not detected with e-MERLIN at 6.9 GHz 
to a 3ct limit of 0.12 mjy/beam. Comparing the 3-a limit 
at 6.9 GHz with the fiux detected by the e-EVN at 1.4 GHz 
gives a limit on the spectral index of a < —1.4, showing that, 
in this case, the spectral index remains steep at higher fre- 
quencies. A search of the catalogues in the NASA/IPAC Ex- 
tragalactic Database (NEDjf] shows that J072212-I-291042 



^ See http: / /www. e-merli n.ac.uk/tech7] for the current capabili- 
ties of the e-MERLIN array. 
^ http://ned.ipac.caltech.edu/ 
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Table 1. Compact sources with a spectral index 031^^^;, < —1.4 observed with the e-EVN. Names correspond to the NVSS catalogue 
names. The coordinates are from the VLA FIRST Survey. The spectral index is calculated from the VLSS and WENSS flux densities, 
and the flux density ratio is calculated from the NVSS and FIRST catalogues. The VLSS, WENSS, NVSS and FIRST flux densities are 
given in mjy. Positions are given in J2000 coordinates. Note that both the WENSS flux density and spectral index of J151229+470245 
were revised after the observations due to an error discovered in one of the catalogues (see text). 



Name 


RA 


Dec 


_,74MHz 
"325 MHz 


oNVSS 
"^1400 MHz 

1400 MHz 


oVLSS 
■^74 MHz 


cWENSS 
■^325 MHz 


oNVSS 
'-^1400 MHz 


cFIRST 
■^1400 MHz 


J072212+291042 


07:22:12.603 


+29:10:41.51 


-1.45 


1.1 


650 


76 


11 


10 


J082916+383453 


08:29:17.351 


+38:34:52.62 


-1.44 


0.9 


940 


112 


11 


12 


J130612+514407 


13:06:12.176 


+51:44:06.93 


-1.58 


1.1 


1640 


159 


27 


25 


J131655+483200 


13:16:55.642 


+48:32:00.13 


-1.52 


1.1 


2390 


252 


67 


61 


J151229+470245 


15:12:29.174 


+47:02:45.59 


-0.50 


1.0 


1130 


543 


217 


221 



has no known counterparts in other wavebands. The com- 
pact nature of this object suggests that the most likely ex- 
planation is either a steep spectrum quasar core or a Galac- 
tic pulsar. Although this source has a low galactic latitude 
(see Table [2]) a search of pulsar catalogues around this lo- 
cation shows no currently known pulsar or pulsar candidate 
at these coordinates (Eatough, private communication). 



3.2 J130612+514407 

Of the detected sources, J130612+514407 has the steepest 
spectral index at low frequencies. Without an observation 
at 6.8 GHz, however, it is impossible to tell if the spectrum 
remains steep at higher frequencies - the positional offset 
in the e-EVN observations was only discovered after the e- 
MERLIN observations had taken place. The e-EVN data 
show an unresolved point source with an integrated flux den- 
sity of 7.1 mJy, recovering ~ 30% of the emission detected by 
the Westerbork array during the same observation, imply- 
ing structure on intermediate scales at 1.4 GHz. However, 
the primary beam of the Westerbork tied-array is signifi- 
cantly smaller than the 12.8" offset of the source position 
from the pointing centre of the e-EVN observation; imaging 
the source without including the Westerbork array results in 
a peak flux density of 18.5 mJy and an integrated intensity 
of 21.6 mJy, recovering almost all of the WSRT-only flux, 
although it should be noted that these measurements are 
somewhat uncertain due to the large positional offset. 

Unlike J072212+291042, J130612+514407 has an opti- 
cal counterpart in SDSS, J130612.15+514407.0, a diffuse- 
looking galaxy with an SDSS spectrum giving a red- 
shift of 0.2773. The source also has a counterpart in the 
JVAS/CLASS sample with an integrated flux density at 
8.4 GHz of 16.2 mJy. If the source does not vary, then this 
results in a spectral index of —0.31 between 1.4 (NVSS) and 
8.4 GHz making it slightly steep, but certainly not ultra- 
steep at these frequencies. However, the CLASS observations 
were carried out with the VLA in its largest A configuration 
resulting in an angular resolution of 220 milli-arcseconds 
(mas), very different to the beam size of the NVSS survey, 
making this spectral index measurement unreliable. This ob- 
ject also hcis counterparts in 2MASS and the ROSAT X-ray 
catalogue, and is cla ssified as the brightest cluster galaxy in 
iKoester et al] |20o3), making the quasar core scenario the 
most likely. 



3.3 J151229+470245 

In contrast, J151229+470245 is clearly resolved in both the 
e-EVN and e-MERLIN images (see Fig. [TJ. The WSRT- 
only data show an unresolved source, as expected from the 
NVSS/FIRST images, with a peak fiux of 192mJy/beam 
and an integrated flux density of 203 mJy. This is within 8% 
of the flux recorded in the FIRST and NVSS catalogues. 
The e-EVN observation shows a source elongated in a NE- 
SW direction, with a peak brightness to total fiux density 
ratio of 1:4. 

The map obtained with e-MERLIN at 6.9 GHz (Fig. 
[1} shows a very different structure. The source is clearly 
resolved and appears to show a core-jet morphology. The 
e-EVN component is coincident with the brightest feature 
observed with e-MERLIN. A two-component fit was carried 
out on the e-MERLIN map using the AIPS task jmfit: a rel- 
atively compact (compared to the e-MERLIN beam), bright 
Gaussian with a peak of 39.5 mjy/beam, an integrated fiux 
of 66.1 mJy and a size of 98x73 mas at a position angle 
of 71 degrees, and a more extended Gaussian "jet" com- 
ponent with a fitted peak brightness of 4.14 mJy/beam, an 
integrated flux density of 28.3 mJy with a size of 0.2x0.12 
arcseconds at a position angle of 145 degrees. 

The nature of this source is rather puzzling. Interest- 
ingly, while the e-MERLIN map (taken on its own) appears 
to show a classic core-jet morphology, the e-EVN ("core") 
component is very resolved suggesting that it is not a sim- 
ple AGN core. Unless there is very strong scatter broaden- 
ing of the target, the observed structure is more consistent 
with lobe emission detected on 10-mas scales. The observed 
steep spectrum, low brightness temperature and no variabil- 
ity (between VLA and WSRT epochs) all support this sce- 
nario. Therefore we may have identifled a "naked jet-lobe" 
source without strong AGN core emission. The absence of 
a compact, self-absorbed core is a sign of a Type 2 AGN, 
when the jet is viewed at a large angle compared to our line 
of sight. It is intriguing though, why such an unbeamed ob- 
ject would have a prominent one-sided jet morphology. An 
alternative interpretation is that we are observing a two- 
sided jet-lobe with a very weak core within the structure. 

On inspection of the WENSS map, we found that 
J151229+470245 is a complex source at the resolution of 
Westerbork at 325 MHz. An examination of the map in com- 
parison with the survey database, the corresponding flelds 
in both NVSS and FIRST, and the new Westerbork data 
obtained during the e-EVN observations, determined that 
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Blue contours: 1512+470 IPOL 688S.675 MHz e-MERLIN 
Red contours: 1512+470 IPOL 1664.490 MHz e-EVN 




15 12 29.21 29.20 29.19 29.18 

RIGHT ASCENSION (J2000) 

Levs = 2.320E-04 * (-1, 1, 2, 4, 8, 16, 32, 64, 128) 
Levs - 3.607E-04 •(-1,1,2, 4, 8, 16) 



Figure 1. e-MERLIN map of J151229+470245 at 6.9 GHz 
(blue contours) together with the e-EVN map at 1.4 GHz (red 
contours). Contours are plotted at -1,1,2,4,8,16,32,64,128 X 
0.23mJy/beam for the e-MERLIN data and -1,1,2,4,8,16 X 
0.36mjy/beam for the e-EVN data. The boxes in the lower left 
indicate the size of the restoring beam for each image. 



the fluxes listed in the WENSS catalogue do not refer to the 
source components at the listed coordinates. The calibrated 
fits file for the WENSS field was retrieved and Gaussian 
components were fitted to the three sources in the region. 
This showed that the source in which we are interested actu- 
ally had an integrated flux of 543 mjy, not 123 mjy as listed 
in the WENSS catalogue. This gives the source a spectral in- 
dex of -0.50 between 74 and 325 MHz, making it steep spec- 
trum, but not ultra-steep. The spectrum is slightly steeper 
between WENSS and FIRST/NVSS with a = -0.62 between 
325 MHz and 1.4 GHz. Therefore, based on its spectrum and 
the observed VLBI and e-MERLIN radio structure we may 
classify J151229-I-470245 as a peculiar compact steep spec- 
trum radio source (CSS). 

This source is also an example of the class of radio 
sources known as Infrared Faint Radio Sources (IFRS), those 
with clear r adio emission but no corresponding infrared de- 
tection (e.g. lNorris et al1l2006l ). Such sources are difficult to 
identify due to their lack of counterparts outside the radio 
regime and, as yet, no IFRS source has a measured red- 
shift, but they are thought to be high redshift quasars due 
to their generally compact nature, mostly steep spectral in- 
dexes and high brightness temperatures. While some of the 
known IFRS sources are extended on arcsecond scales with 
the VLA, and a small number of such source s have previ- 
ously been detected at VLB I-resolution (e.g. iNorris et ahl 
120071 : iMiddelberg et al. I [20081 ). this appears to be the only 
example so far which is resolved on VLBI scales. 



4 DISCUSSION AND CONCLUSIONS 

Of the five sources in the sample, one, J151229-f470245, is 
not an ultra steep spectrum source. Of the remaining four. 



two are undetected with VLBI, two are detected but are un- 
resolved on milli-arcsecond scales. These new high resolution 
data are not sufficient to uniquely determine the nature of 
these sources, however they do provide valuable information 
on the variety of sources one may find in a larger USS sam- 
ple. As one may expect, a significant fraction of these sources 
are completely resolved out since steep spectrum emission 
is often related to very extended lobes in radio galaxies. A 
detection rate of fifty per cent, however, indicates that very 
high resolution studies of a moderate fraction of USS sources 
will be possible. 

Both of our USS source detections, J072212-)-291042 
and J130612-I-514407, are consistent with a point source 
on VLBI scales. Note that in these short, exploratory ob- 
servations with limited uu-coverage (especially on the long 
baselines) we cannot adequately constrain the source sizes, 
but they do appear more compact than our beam of 32x26 
mas. The emission in these two cases is, therefore, consistent 
with either a Galactic pulsar or a very steep spectrum quasar 
core. The former hypothesis will be easily tested by further 
VLBI observations that can potentially show the proper mo- 
tion of the pulsar in our Galaxy. 

Since there are generally no detections of these sources 
in any other wavelength regime and no redshift information 
is available, they are very difficult to classify. These new 
radio observations provide fresh information and vital clues 
in the search for an understanding of their nature. 

J072212-I-291042 is compact on VLBI scales but the e- 
EVN observations presented here recover less than 20% of 
the flux seen at the same frequency on VLA scales, clearly 
indicating the existence of structure on intermediate spa- 
tial scales to which neither of these arrays are sensitive, e- 
MERLIN observations at 6.9 GHz fail to detect the source to 
a 3-cr limit of 120 ^Jy/beam, confirming the steep spectrum 
extends to higher frequencies. Observations with e-MERLIN 
at 1.4 GHz would probe this intermediate-scale emission, but 
this band was not available at the time when the commis- 
sioning observations described here were obtained. The fact 
that only ~20 per cent of the FIRST/NVSS flux is recov- 
ered on VLBI scales strongly suggests the existence of struc- 
ture on interm e diate scales for these sources. A study by 
iKlockner et al.l (|2009l ) of a sample of 11 high-redshift ob- 
scured quasars also found similar evidence for structures 
on intermediate scales, although the VLBI-recovered fluxes 
there ranged from 30 to 100 per cent of the low-resolution 
measurements. 

J130612-I-514407 is also compact on VLBI scales with 
significantly less flux lost when imaged with long baselines. 
Unlike J072212-I-291042, however, this source does have 
counterparts in other wavebands which suggest the source 
is associated with a galaxy at a redshift of 0.2773. 

While J151229-I-470245 appears to be a compact steep 
spectrum object, rather than an ultra steep spectrum source, 
it is still compact with a steep spectrum at low frequen- 
cies and interesting structure at higher resolutions. VLBI 
observations of J151229-I-470245 only recover 20% of the 
NVSS flux. In sharp contrast to J072212-I-291042 and 
J130612-I-514407 however, the object is clearly resolved on 
both e-MERLIN and VLBI scales, showing different struc- 
tures in the two maps. While the e-MERLIN map shows 
what appears to be a classic core-jet morphology, the "core" 
is clearly resolved in the e-EVN map. Observations with e- 
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Table 2. e-EVN detected sources. The measured coordinates in the VLBI experiment (±0.1 mas) and the corresponding galactic coordi- 
nates are both given. Flux densities are given for each observation: WSRT (1.6 GHz), e-EVN (1.6 GHz) and e-MERLIN (6.8 GHz), and 
the NVSS-WSRT spectral index, and VLBI- WSRT flux density ratios are calculated. For J072212H-291042 a 3a upper limit is given at 
6.89 GHz. Note that the WSRT and VLBI fluxes at 1.6 GHz were obtained simultaneously, and that J130612-I-514407 was not observed 
with e-MERLIN. 



Name 


RA 


Dec 


m 


oWSRT 
"-^leTOMHz 


cc-EVN 
'-^1670MHz 


CO-MERLIN 

■-^esgoMHz 


^1400 MHz 
"1670 MHz 


ge-EVN 
'-'1670MHz 
cWBRT 
1670MHz 


J072212+291042 


07;22;12.60906 


+29:10:41.6648 


189, 19 


7.6 


2.1 


<0.12 


-1.82 


0.28 


J130612+514407 


13;06;12. 16262 


+51:44:06.9627 


118, 65 


23.1 


21.6 




-1.28 


0.94 


J151229-I-470245 


15:12:29.19905 


+47:02:45.3837 


78, 56 


203 


42.0 


94.0 


-0.51 


0.21 



MERLIN at L-band (1.4 GHz) would again allow us to probe 
the morphology of the missing flux. 

While firm conclusions regarding the nature of these 
sources remain elusive, the observations and unexpected re- 
sults presented here do show that comijact ultra-stccp spec- 
trum sources are a diverse class of object whose nature can 
be probed with VLBI. 
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